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The stability of the equilibrium configurations of a nematic liquid 
crystal confined between two coaxial cylinders is analysed when a 
radial electric field is applied and the flexoelectric effect is taken 
into account. The threshold for perturbations depending only on the 
radius r in the cylindrical coordinate system and strong boundary 
conditions is studied. A new type of orientational transition caused 
by pure flexoelectric effect is found. 
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INTRODUCTION 

The uniform orientation of nematic liquid crystal (NLC) is caused by 
the orienting action of appropriately treated confining plates, which 
define a fixed orientation of the NLC at the boundary. When an elec- 
tric or magnetic field is applied to such NLC sample in an appropriate 
direction, a director reorientation via a Freedericksz transition may 
result, which has been studied extensively in case of a plane layer 
geometry. 

For a NLC confined between two coaxial cylinders the orienta- 
tional transition may occur even in the absence of external fields 
depending on the ratio of the radii of the inner cylinder to that of 
the outer one |2[. The Freedericksz transition in this geometry 
was analysed for a magnetic field directed radially outwards from 
the common axis of the cylinders and tangentially to circles in 
the plane perpendicular to the cylinders' axis @, |J as well as for 
a radial electric field. Here a voltage was applied between the two 
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coaxial cylinders |||, [?]]. In this geometry the critical voltage for the 
Freedericksz transition depends on the ratio between the cylinders 
radii. 

In this paper we analyse the influence of the flexoelectric effect on 
the orientational transitions in a NLC confined between two coaxial 
cylinders under applied d.c. radial electric field. The flexoelectricity 
describes the linear coupling between an applied electric field and 
gradients in the director field ||. In contrast to the case of a plane 
nematic layer, in cylindrical geometry under applied electric field one 
has a contribution from the flexoelectric effect in the bulk torque act- 
ing on the director. Here we calculate the critical voltage for the ori- 
entational transition where the director distribution n = (n r ,n^,,n z ) 
depends only on the radius r in the cylindrical coordinate system 
(r, </>, z) with z axis along the axis of cylinders. Strong boundary 
conditions are imposed. Two typical initial director distributions 
[homeotropic orientation n = (1,0,0) and planar one n = (0,1,0)] 
are considered. It is found that for planar orientation a new type of 
orientational transition caused purely by the flexoelectric effect takes 
place. 



BASIC EQUATIONS 

Let us consider the nematic enclosed between two infinite coaxial 
cylindrical electrodes whose inner and outer radii are r\ and Ti, re- 
spectively. The director distribution n = n(r) depends only on r. 
The voltage applied between the electrodes provides the radially di- 
rected electric field E = (E,0,0). The equilibrium director configu- 
ration can be found from the extremum condition of the total free 
energy (per unit length of cylinders) 

T = 2ix \ (F d + F el + F fl )rdr . (1) 

Here F d is the elastic energy density 

F d =- 

The electric field contribution is given by 



K x (V • n) 2 + K 2 (n • (V x n)) 2 + K 3 (n x (V x n)) 2 . (2) 



F el = ~e e a (h-E) 2 , (3) 
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where e a = en — ej_ is the dielectric anisotropy and the radial compo- 
nent of the electric field E = U/[r ln(r 2 /r!)] with U = U(r{) — Uir-i) 
the voltage drop across the cylinders. The contribution of flexoelec- 
tricity 

F fl = - (E ■ P) , P = e n n (V • n) - e 33 (n x (V x n)) (4) 

is connected with flexo-polarization P induced by the splay and bend 
deformations of the director field and en, 633 are the flexoelectric 
coefficients. 

Minimisation of (JTJ) with the normalisation n 2 = 1 gives the 
torque balance equation 

[fixh] = ' * = i{£$-^> i=T >*>* (5) 

where / = (F d + F ei + Ffi)r. One can easily verify that n = (1, 0, 0) 
and n = (0, 1,0) are solutions of equation (|) for homeotropic and 
planar boundary conditions, respectively. 



STABILITY ANALYSIS 



Let us now examine the stability of these configurations. 
I) Homeotropic anchoring with boundary conditions 

(1,0,0). 



nir 



nir 



r 2 



Small deviations from the exact solution no = (1, 0, 0) of 
written as 

1 



n = n + Sh 



-(Sn^ + 8n z ); Sn^; 8n 2 



(6) 
can be 

(7) 



up to second order in the perturbations Sn^, 5n z which depend only 
on r. Then for the change of total free energy one obtain up to 
second order in 5h 



= 27r £ 2 {^ l[ -^K + 8*1)} + \ K ^< + 5n * + ^ 
Ue n — (8nl + 5n 2 z ) + Ue m —5n 2 Ardr , 



(8) 
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where the boundary conditions 5h(r — 7*1) = 5h(r = r 2 ) = were 
taken into account, U = C//hi(r 2 /ri) and /' = df /dr. 

The necessary and sufficient condition for AjF feom to be positive 
definite is that the Euler-Lagrange equations 

d 2 ~ 
K3 d^ Sn * + [Kl ~ Ks ~ tQtaU2 ~ 2C/(en + e ^ 5n t> = ' (9) 

K 3^n z + [K, - e e a U 2 - 2Ue n ]5n z = , (10) 

have nonzero solutions 5ns, 5n z 7^ in the range < a < ln(r 2 /ri) 
where we introduced the new variable a = ln(r/ri). Since equations 
(0), ([TDD are uncoupled one gets two independent conditions for the 
instability of initial homeotropic orientation with respect to 5ns and 
5n z perturbations. Solving (Q), flipp one obtains for the critical volt- 
age U c for the two types of orientational transitions in cylindrical 
geometry 



5ns : -sign(e a 



U C V ( U c \ 2(e 11 + e 33 )ln(r 2 /r 1 ) 



U F J \U F J nJK 3 e \e a 



K 3 -K x f\n(r 2 / ri )\ , s 



5n 2 : -sign(e 



f/ c \ 2 /C/ c \ 2e 11 ln(r 2 /r 1 ) 



f/ F3 y V^W 7r,/^ 3 eo|e a 



where Up 3 = ir^J K 3 /e \e a \ is the threshold voltage for the bend Free- 
dericksz transition in case of a plane NLC layer. In Figure [I] the 
two types of transitions correspond to 5ns and 5n z perturbations are 
shown schematically. In the absence of electric field (U = 0) there 
is a possibility of orientational transition for a radius ratio above 
some critical value. The transition corresponding to 5n z perturbation 



[Fig. 1(b)] comes first and occurs at ln(r 2 /ri) = n y K 3 / K\. 

For a nematic liquid crystal with negative dielectric anisotropy 
(e a < 0) the flexoelectric effect reduces the critical voltage for the 
orientational transition. Using the MBBA material parameters at 
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Figure 1: Orientational transition for 611^ (a) and 5n z (b). 




Figure 2: Critical voltage versus ratio of the radii for MBBA with 
homeotropic (a) and planar (b) anchoring. 



25°C K x = 6.66 ■ 1(T 12 N, K 2 = 4.2 ■ KT 12 N, K 3 = 8.61 • 1(T 12 N, 
e a = —0.53 and taking for the flexocoefficients the order of magnitude 
-10" 11 C/m " 



10 we found that the transition 



en + e 33 = e n -- 

corresponding to 5n z perturbation has lower threshold [Fig. @(a)]. 
In this case above r 2 /ri pa 35.6 the orientational transition occurs 
without electric field. 



II) Planar anchoring with boundary conditions 
n(r = 7"i) = n(r = r 2 ) = (0, 1, 0) 



(13) 



Similar to the case of homeotropic orientation, small deviations from 
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Figure 3: Orientational transition for 5n r (a) and 5n z (b). 



the exact solution no = (0, 1, 0) of @ can be written as 

1 



n = n + dn 



Sn r ; 1 — -(Sn 2 + 5n 2 z ); 5n^ 



(14) 



up to second order in the perturbations 5n r , 5n z which depend only 
on r. Then the change of total free energy up to the second order in 
5h is 

kF p i an = J^(n ) - ^(no + 8h) = 

pj"2 "I 1 1 1 

= 2vr / {-^[K 2 + -M\ + 9^[<K 2 + -M] - 

J r\ Z T Z T 

11 1-1 

~K 3 — (8n 2 r + 5n 2 z ) - -e e a U 2 —5n 2 r - 
Z r z r 

1 1 

Uen—Sn 2 . - Ue 33 — (5n 2 + 5n 2 z )}rdr , (15) 

where the boundary conditions 5h(r = ri) = 5h(r = r 2 ) = were 
taken into account. 

The instability conditions can be obtained from the Euler-Lagrange 
equations 

d 2 ~ 
Kx—Jnr + [K 3 -K x + e e a U 2 + 2U(e u + e 33 )}5n r = , (16) 



da 2 
d 2 

'da 



K 2 —5n z + [2K 3 -K 2 + 2Ue 33 }Sn z = . 



(17) 
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Here again the equations ([Rj)> fllTD are uncoupled and one has two 
independent conditions for the instability of initial planar orientation 
with respect to 5n r and 8n z perturbations (Fig. |3|). The critical 
voltage U c for the orientational transition can be found from the 
solution of (|T|), flT7|) 

U c \\ ( U c \ 2(eu + e33)ln(r 2 /ri) 



6n r : sign(e a ) — + 



fJ 7rVKie |e a 

2 



K 3 -K, An(r 2 /n) \ 

K X \ TT ) 



(18) 



U c \ 2e 33 ln(r 2 /ri) 



F 2 y 7r,/K 2 eo|e a | 



2^3-^2 / ln^/rQ V 
i^2 V ^ ) 



(19) 



where Up 1 = rcy Ki/e \e a \, Up 2 = K 2 /e \e a \ are the threshold 
voltages for the splay and twist Freedericksz transition in the case 
of a plane NLC layer, respectively. Without electric field one has 
the orientational transitions at the critical radius ratio ln(r 2 /r!) = 



n^K./iKs - K,) [Fig. |(a)] and ln(r 2 /n) = 7r^/K 2 /(2K 3 - K 2 ) 
[Fig. |(b)]. 

For NLC with e a > the flexoelectric effect reduces the critical 
voltage for the transition corresponding to Sn r perturbations similar 
to that found in the case of homeotropically oriented nematics with 
e a < [Fig. 0(a)]. Equation fllTf ) demonstrates a new feature of 
the cylindrical geometry, namely, the possibility of an orientational 
transition caused purely by the flexoelectric effect. In Figure ^|(b) 
the threshold voltage corresponding to the flexoelectric instability 
(Sn z perturbation) as a function of the ratio of the radii is plotted 
for the MBBA material parameters and = — 10~ n C/m. Above 
r 2 /ri w 5.96 the transition takes place without electric field. 

Thus, the analysis of the orientational transitions in a NLC con- 
fined between two coaxial cylinders shows the strong influence of the 
flexoelectric effect on the critical voltage. For planar boundary condi- 
tions a pure flexoelectric instability is found. In this case the polarity 
of critical voltage depends only on the sign of the flexo-coefficient 633 
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which allows to determine the value and the sign of e 33 in a simple 
experiment. Choosing the radius ratio of the cylinders close to the 
critical (above which the transition occurs without electric field) one 
can make the threshold voltage for the flexoelectric-induced orienta- 
tional transition very small and reduce the possible influence of the 
electric current in a not sufficiently clean NLC materials. 
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